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SUMMARY 


INTRODUCTION 

One  of  the  objectives  of  TED  Project  No.  ADC  AV-4401G,  reference 
(a),  was  to  provide  cooling  specifications  for  airborne  electronic  equip¬ 
ment.  Report  No.  NADC- EL-59113,  reference  (b),  delineated  the  tech¬ 
nical  requirements  upon  which  the  specification  should  be  developed. 

It  recommended  an  analytical  specification,  the  requirements  of  which 
should  be  based  upon  a  unit  operations  philosophy. 

The  proposed  program  for  the  development  of  an  analytical  cooling 
specification  was  accepted  by  the  Bureau  of  Naval  Weapons  (BUWEPS), 
reference  (c),  as  being  satisfactory  and  recommended  that  it  be  carried 
out  as  submitted  in  reference  (b). 

The  technical  elements  of  the  proposed  analytical  specification  will 
call  out  the  rules  that  require  the  contractor  to  analytically  prove  that 
his  cooling  package  design  will  provide  reliable  component  temperature 
levels  within  the  electrical  equipment  and  maintain  cooling  hardware 
penalties  such  as  weight  and  fuel  consumed  in  overcoming  momentum 
drag  of  the  aircraft  to  a  minimum.  Reference  (b)  pointed  out  that  all 
participants  involved  either  in  the  preparation  or  review  of  the  reports 
would  require  analytical  design  tools  for  determining  whether  tempera¬ 
ture  and  penalty  values  are  being  kept  to  a  minimum.  Analytical  tools 
of  heat  transfer  for  determining  the  performance  characteristics  of 
various  types  of  unit  operations  of  cooling  as  well  as  weight  penalty 
factors  are  to  be  covered  in  a  series  of  reports. 

The  first  report  on  analytical  tools  deals  with  the  unit  operation  of 
cooling  known  as  regenerative  heat  sinks  for  airborne  electronic 
equipment. 

Regenerative  heat  sinks  are  a  means  of  removing  heat  from  a  heat 
sink  in  electronic  equipment  during  a  low  speed  phase  of  a  flight  mis¬ 
sion.  During  the  high  speed  phase  of  the  mission,  the  electronic  equip¬ 
ment  can  maintain  component  temperatures  at  reliable  levels  due  to  the 
Initial  low  temperature  and  heat  capacity  characteristics  of  the 
regenerator. 
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SUMMARY  OF  RESULTS 

The  effectiveness  of  an  analytical  specification  will  be  dependent 
upon  the  availability  and  utilization  of  analytical  tools  of  design  covering 
various  unit  operations  of  cooling  airborne  electronic  equipment  and 
methods  of  calculating  penalty  factors  to  the  over-all  aircraft  system. 

The  packed  bed  regenerator  is  considered  to  be  an  Important  unit 
operation  of  equipment  cooling  design.  Analytical  theory  and  equations 
have  been  provided  which  can  be  used  for  the  design  of  packed  bed 
regenerators.  Additional  analyses  required  for  different  and  more  com¬ 
plicated  types  of  boundary  conditions  will  be  covered  in  subsequent 
reports. 

CONCLUSIONS  AND  RECOMMENDATIONS 

Methods  of  analysis  have  been  provided  which  allow  calculations 
to  be  made  for  the  size  and  weight  of  packed  bed  regenerator  required 
to  maintain  temperature  limits  within  a  required  range.  Future  reports 
will  cover  unalytical  tools  for  converting  weight  of  packed  bed  regenera¬ 
tors  to  over-all  penalty  factors  for  aircraft  systems. 

It  is  recommended  that  this  report  be  made  available  to  navy  con¬ 
tractors  for  information  purposes  and  to  be  used  as  an  aid  in  their 
technical  development  programs. 


( 
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DISCUSSION 


GENERAL 

Electronic  engineers  for  some  time  have  been  depending  on  the 
principle  of  regeneration,  known  better  as  ’’thermal  inertia,"  for  main¬ 
taining  components  of  equipment  packages  at  reliable  levels  during 
supersonic  flight.  Generally  for  certain  aircraft  systems,  it  was  rec¬ 
ognized  that  no  cooling  or  an  inadequate  amount  of  refrigerated  cooling 
air  would  be  supplied  to  the  electronic  equipment  during  supersonic 
flight.  The  heat  capacity  of  the  mass  of  the  electronic  equipment  had 
to  be  depended  upon  to  absorb  heat  from  electronic  components  for  a 
given  period  of  time. 

In  general,  no  attempt  was  made  to  calculate,  prior  to  hardware 
development,  temperatures  of  the  equipment  or  whether  penalty  factors 
would  be  maintained  at  a  minimum.  Further,  since  analytical  tools  of 
design  did  not  exist,  it  was  not  possible  to  analytically  determine 
whether  a  regenerative  heat  sink  or  another  type  of  cooling  should  have 
been  used  for  a  particular  operation  of  a  flight  mission. 

The  discussion  presented  below  provides  theory  and  mathematical 
equations  for  calculating  temperatures  for  cooling  and  heating  packed 
bed  regenerators. 

THEORY  OF  REGENERATION  OF  PACKED  BEDS 

The  temperatures  of  the  air  and  packed  bed  of  a  regenerative  heat 
exchanger  depend  on  space  and  time.  Only  certain  analytical  proce¬ 
dures  are  considered  in  this  paper.  Some  of  the  main  publications  of 
analytical  solutions  are  based  on  the  studies  of  Heiligenstaedt,  Hausen, 
Nusselt,  and  Schumann,  references  (d),  (e),  (f),  and  (g).  One  of  the 
most  important  of  the  analytical  solutions  is  attributed  to  Schumann 
who  solved  the  complicated  theoretical  heat- transfer  rate  equations  for 
the  simple  case  of  an  incompressible  fluid  passing  through  a  bed  of 
solid  particles.  The  following  equation 

-  (dT/dY)  =  T-0  =(d9/d  Z) 
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describes  the  heating  or  cooling  of  a  packed  bed  with  the  initial  condl 
tions 


Y  =  0; 

Z  =  0; 


T  =  T 
0  =  0 


by  the  passage  of  a  fluid  entering  at  a  temperature  T0.  T  and  9  are 
the  temperatures  of  the  fluid  and  packed  bed,  respectively,  Y  units 
from  the  inlet  plane  at  a  time  given  by  the  group  Z,  where 

Y  =  (hX/>/CG),  dimensionless 

h  =  heat-transfer  coefficient,  BTU/hr  ft3  °F 

X  =  distance  from  entrance  of  bed,  ft 

P  =  fluid  density,  lb/ft3 

C  =  fluid  heat  capacity  per  unit  volume  at  constant  pressure, 
BTU/ft3  °F 

G  =  fluid  flow  rate  based  on  a  superficial  cross-sectional  area 
of  packed  bed,  lb/hr  ft3 

Z  =  {h£t  -  (X/v)jJ  ,/Cb(l  -  f),  dimensionless 

t  =  time  after  start  of  fluid  flow,  hr 

v  =  fluid  flow  per  superficial  cross-sectional  area,  ft3/ft3  hr 
0^  =  packed-bed  heat  capacity  per  unit  volume,  BTU/ft3  °F 
f  =  void  fraction,  dimensionless 


In  his  analysis,  Schumann  assumed: 

1.  The  particles  were  so  small  or  had  such  a  high  thermal  dif- 
fusivity  that  any  given  particle  may  be  considered  to  be  at  a  uniform 
temperature. 


( 
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2.  The  resistance  to  transfer  of  heat  by  conduction  in  the  solid 
itself  was  negligible. 

3.  The  rate  of  heat  transfer  from  the  fluid  to  a  solid  at  any  point 
was  proportional  to  the  average  difference  in  temperature  between 
fluid  and  solid  at  any  given  point. 

4.  The  change  in  volume  of  the  fluid  and  solid  due  to  the  change 
in  temperature  was  negligible. 

5.  The  thermal  constants  were  Independent  of  temperature.  With 
the  Initial  temperature  of  the  bed  equal  to  zero,  he  produced  the  fol¬ 
lowing  results: 

(0/To)  =  1  -  e-y-z  £  YnMn(YZ)  =  e‘y"z 

n  =  0 

(T/T0)  =  1  -  e'y~z  £  YnMn(YZ)  =  e  'y"z 

n  =  1 

where  Mn(YZ)  is  the  nth  derivative, 

[dnM0(YZ)]  /  [d(YZ)n] 

with  respect  to  YZ,  and 

M0(YZ)  =  1  +  YZ  +  [  (YZ)2/2 1]  +  [  (YZ)3/3 !  ]  +  .  .  .  . 

Schumann  provided  numerical  solutions  in  the  range  Y  =  0  -  10, 

Z  =  0  -  9|  in  the  form  of  graphs. 

DEVELOPMENT  OF  EQUATIONS  FOR  COOLING  OF  PACKED  BEDS 

There  appears  to  be  little  published  literature  for  using  numerical 
solutions  with  a  digital  computer  to  determine  the  time  required  to  cool 
the  packed  bed  with  power  being  dissipated  as  heat  into  the  packed  bed. 
In  the  development  of  the  method  of  finite  differences,  for  cooling  the 
packed  bed,  the  annular  cylinder  containing  the  bed  material  is  divided 
into  a  number  of  equally  spaced  intervals  each  of  lengthAS,  as  shown 
in  figure  1. 


f  ZnM  (YZ) 
n  =  1  n 

f  Z"Mn(YZ) 
n  -  0 
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An  energy  balance  is  made  on  an  elemental  volume.  This  total 
energy  includes  a  term  Hs  for  rate  of  external  heat  addition  (by  the 
electronic  equipment)  to  the  packed  bed  per  unit  volume.  If  no  heat  is 
being  added  to  the  packed  bed,  as  might  occur  with  the  electronic  equip* 
ment  in  a  standby  position,  the  heat  term  H„  is  dropped  from  equation 

a). 


HgAdS  +  p  Cp(dT/dt)AdS  = 

GoCp(dT/0S)AdS  +  (1  -  f),bCb(d0/dt)AdS  (1) 

where 

Cp  =  specific  heat  of  air,  BTU/lb  °F 

Af  =  cross  sectional  area  of  bed,  ft^ 

dS  -  incremental  height  of  bed,  ft 

Pfr  =  specific  heat  of  packed  bed,  lb/fr 

G0  =  flow  rate  per  unit  cross  sectional  area,  lb/hr  ft^ 

Since  pCp(dT/<3t)AdS  is  considered  negligible  compared  to  HgAdS  or 
G0Cp(  T/d  s)AdS,  equation  (1)  becomes 

HsAdS  -  GoCp(0T/0S)AdS  =  (1  -  f)  p  ^9  /d  t)AdS  (2) 

Now  it  is  to  be  noted  that 

(1  -  f)PbCb(00/dt)AdS  =  ha(T  -0  )AdS 

and  equation  (2)  becomes 

Hs  -  G0Cp(aT/dS)  =  ha(T  -0  )  (3) 

To  transpose  equation  (3)  in  terms  of  finite  differences,  let 
dT/aS  =  (Tra+i  n-Tmi  „)/AS; 

dfl/at  =  (em+i,  „-em> 

^m+(l/2),  n  “  ^m+1,  n  +  ^m, 
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Now  substitution  of  the  finite  difference  terms  into  equation  (3)  results 
in 


Hs  -  G0Cp  <Tm+1>  „  -  Tro>  n)/AS  = 

ha(Tm+ (1/2),  n  "  9  m+(l/2),  n>  <4) 


rearranging  terms, 


H_  GnC  T  ,  n  GnC_T 
s  op  m+1,  n  0  p  m,  n 

ha  haAS  haAS 


(Tm+(l/2),  n  "  B  m+(l/2),  n* 


By  using  a  modulus  M  =  (G0C  /haAS)  and  substituting  into  equation  (5), 
there  results  p 


(6) 


(Ho/ha)  -  MT  .  +  MT  = 

'  s'  '  m+l^  n  m,  n 

^Tm+(l/2),  n~B  m+(l/2),  n) 

again  rearranging  terms  gives 

T  -  T  = 

m+1,  n  m,  n 

(l/M)  [flm+(1/2)>  n  -  Tm+(l/2),  n  +  ,Hs/ha>] 
substituting  Tm+(1/2)i  „  gives 

Tm+l,  „  ■  [»/0  ♦  no]  x  [2®lr+(1/2)i  „  ♦  Tm>  „  (2M  -  1)  ♦  k] 


(7) 


(8) 


where  K  =  2Hs/ha.  It  is  to  be  noted  that  for  the  case  of  no  power  being 
dissipated  as  heat  into  the  packed  bed,  the  term  K  =  2Hs/ha  is  taken 
out  of  equation  (8).  A  finite  energy  balance  is  now  made  on  both  the  air 
and  the  bed  material  contained  in  an  elemental  volume. 
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GoCpTm,  „  '  GocpTm+l,  n 


c&s/At)(i  -  f)pbcbem+(1/2)(  n  +  j  ■ 
-  o,bcbem+(1/2)>  n 


Further  manipulation  gives 
GoCp(Trr,  n  -  Tm+1,  n>  = 


(9) 


(AS/At)(l  -  f)pbCb  ^m+(i/2),  n+l  "  8  m+(l/2),  n)  (10) 

Now  let  a  modulus  N  =  |AS  (1  -  f)pbCb] /(G0CpAt),  equation  (10)  becomes 
^m+(l/2),  n+l  “  ^m+(l/2),  n  =  (l/NHT,^  n  -  Tm+^  n)  (11) 


Equations  (8)  and  (11)  are  now  combined.  For  the  case  of  no  heat  added 
to  the  bed,  this  becomes 


^m+(l/2),  n+l  =  (1/2) {[l  -  2M  + 

(2/N]Tm>  n  +  [  1  +  2M  -  (2/N)]Tm+i  J.  (12) 


Equations  (8)  and  (12)  become  the  working  equations  for  determining 
the  temperature  of  a  packing  material  and  air  with  respect  to  time  and 
position  within  the  packed  bed. 


It  is  necessary  to  select  suitable  values  of  M  and  N  consistent  with 
the  physical  configuration  under  investigation.  Examination  of  equation 
(8)  for  M  >  (1/2)  discloses  that  a  negative  effect  on  Tm  n  is  avoided  if 
M  =  1.  Suitable  N  values  may  be  estimated  if  Tm+j  is  efiminated 
between  equations  (8)  and  (12).  The  value  selected  for  M  fixes  A  S;  the 
value  selected  for  N  then  fixes  At  when  AS  is  fixed. 


For  example,  when  M  and  N  =  1,  then  equations  (8)  and  (12)  become 
"^m+l,  n  =  (l/2)(20jr+(l/2),  n  +  Tm,  n^  (12) 
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®m+(l/2),  n+1  =  U/2)(Tm>  n  + 

A  similar  approach  is  used  to  cover  the  condition  of  power  dissipated 
as  heat  going  into  the  packed  bed.  Equation  (13)  becomes 

^m+1,  n  ~  U/3)(2®Ir+(i/2)>  n  +  ^m,  n  +  ^)  (^) 

DEVELOPMENT  OF  EQUATIONS  FOR  HEATING  OF  PACKED  BED 

A  discussion  now  follows  on  the  packed  bed  of  the  electronic  equip 
ment  previously  cooled  to  a  predetermined  low-temperature  value  and 
now  acting  as  an  absorber  of  heat  from  the  electronic  components  (the 
entire  equipment  being  sealed  and  insulated  from  a  high-temperature 
compartment  environment  caused  by  high  Mach  Number  flight  of  the 
aircraft). 

The  interstices  between  the  packings  are  filled  with  air.  The  bed 
and  air  for  the  calculations  performed  are  initially  at  a  uniform  tem¬ 
perature.  Heat  from  the  electronic  components  is  conducted  into  the 
container  through  the  walls  as  indicated  in  figure  1.  The  container 
wall  temperature  is  specified  in  terms  of  a  given  temperature-time 
history.  It  is  also  assumed,  for  any  time,  that  the  packing  and  air 
temperatures  at  any  axial  coordinate  are  constant  over  the  cross  sec¬ 
tion  and,  further,  that  the  effects  of  axial  conduction  and  free  convec¬ 
tion  of  are  minor  and  may  be  neglected.  A  heat  balance  on  the 
system  is 

hwgVew  •  T)  =  />Cpf(dT/dt)  +  pbCb(l  -  f)(dS/dt)  (16) 

where 

hWg  =  film  coefficient  between  wall  and  air,  BTU/hr  ft^  °F 

Aw  =  surface  area  of  wall,  ft^ 

0W  =  temperature  of  wall,  °F 

Using  the  same  assumptions,  a  heat  balance  may  be  set  up  for  the  flow 
of  heat  from  air  to  bed  resulting  in  the  following  equation: 
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ha(T-0)=,bCb(l  -f)(d0/dt)  (17) 

Equation  (16)  then  becomes 

hwgAw(0w  “  T)  =  P  Cpf(dT/dt)  +  ha(T  -  9  )  (18) 

Equations  (17)  and  (18)  may  be  conveniently  solved  by  means  of  an  ana¬ 
log  computer. 

Computer  Program 


With  the  equations  derived  above,  the  program  was  divided  into 
three  phases  specified  below: 

1.  In  phase  I,  for  a  given  airflow  rate  and  low  air  temperature 
occurring  during  the  low- speed  phase  of  the  flight  profile,  calculations 
were  made  to  determine  the  time  required  to  lower  the  temperature  of 
the  packed  bed  to  the  temperature  of  the  incoming  air,  with  the  equip¬ 
ment  not  operating. 

2.  Phase  II  was  essentially  the  same  as  phase  I  with  the  exception 
ihat  various  amounts  of  power  dissipated  by  the  electronic  package 
would  be  transferred  into  the  area  of  the  packing. 

3.  In  phase  III,  the  entire  electronic  package  was  considered 
thermally  insulated  from  the  aircraft  compartment  with  no  air  entering 
or  leaving  the  equipment.  All  of  the  heat  dissipated  by  the  electronic 
compartments  was  being  taken  by  the  packings.  The  entire  electronics 
package  was  viewed  as  being  thermally  insulated  from  the  high- tem¬ 
perature  aircraft  environment  by  means  of  insulating  barriers.  The 
packing  materials  were  all  considered  to  be  spherical  with  a  diameter 
of  1/8  inch.  The  materials  were  selected  to  cover  a  range  of  proper¬ 
ties  as  shown  in  table  I. 

For  phase  III,  the  interstices  between  the  packings  are  filled  with 
air.  The  bed  and  air  for  the  calculations  performed  are  initially  at  a 
uniform  temperature.  Heat  is  introduced  into  the  container  wall.  This 
is  sp  cified  in  terms  of  a  given  temperature-time  history. 
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The  physical  configurations  of  the  various  packed  beds  are  shown 
in  table  II.  (All  dimensions  are  in  inches.) 

The  cooling  air  conditions  are  as  specified  in  table  HI.  Two  sets 
of  calculations  were  made,  one  with  an  initial  packed-bed  temperature 
of  86°  F,  and  the  other  140°  F. 

Equations  (13)  and  (14)  were  used  for  phase  I.  A  typical  set  of 
calculations  for  M  =  N  =  1  is  shown  in  figure  2.  Figure  3  shows  curves 
for  phase  II  with  power  dissipated  as  heat  into  the  packed  bed. 

Figure  4  shows  typical  curves  obtained  for  phase  III.  In  examina¬ 
tion  of  the  data  shown  in  figure  4,  it  should  be  realized  that  two  film 
coefficients  of  heat  transfer  are  involved  in  the  problem.  There  is 
the  film  coefficient  of  heat  transfer  between  the  concentric  cylinder 
and  the  air  designated  as  hWgAw  and  between  the  air  and  the  bed  mate¬ 
rial  ha.  During  actual  environmental  flight  conditions  these  may  vary, 
depending  upon  variables  such  as  the  cooling  design  of  the  electronic 
equipment,  the  type  of  aircraft,  or  its  mission. 

Data  were  obtained  to  cover  a  broad  range  of  film  coefficients. 
Figure  4  gives  only  the  extreme  ranges  of  film  coefficients. 


FIGURE  1  -  Regenerative  Packed  Bed  for  Electronic  Equipment  Cooling 
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TABLE  I 

PACKING  MATERIALS  AND  PROPERTIES 


Steel 

Balls 

Glass 

Beads 

Gravel 

Foam 

Density  (lb/ft^) 

497 

157 

165 

10 

Specific  heat  (BTU/lb  °  F) 

0.106 

0.20 

0.65 

0.206 

Surface  area/unit  volume 

334 

334 

334 

334 

Voids 

0.420 

0.420 

0.420 

0.420 

TABLE  II 

PACKED  BED  DIMENSIONS 
(IN) 


Column  Inside  Diameter  of  Outside  Diameter  of 

Length  Concentric  Cylinder  Concentric  Cylinder 


8 

12 

16 

8 

12 

16 

8 

12 

16 


5 

5 

5 

5.5 

5.5 

5.5 

6 
6 
6 
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TABLE  III 

AIR  SUPPLY  CONDITIONS 


Flow  Rate  (lb/hr)  Air  Temperature  (°  F) 


781 

68 

563 

32 

404 

-4 

266 

-40 

175 

-63 

108 

-63 
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FIGURE  2  -  Phase  I  Bed  Temperature- Time  History  for  Foam  Material 
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FIGURE  3  -  Phase 'ultecf^emperature  Time  History 
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FIGURE  4  -  Phase  HI  Bed  Temperature- Time  History  for  Two  Materials 
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